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T—l • ABSTRACT 

> ; 

■ Context. One of the most outstanding problems related to numerical models of galaxy formation is the so-called "angular momentum catastro- 
' phe", i.e., the inability to explain theoretically the high angular momentum observed in local disk galaxies. 
( Aims. We study the evolution of the angular momentum from z~0.6 to z=0 to further our understanding of the mechanisms responsible for 
^ ■ the large angular momenta of disk galaxies observed today. This study is based on a complete sample of 32, 0.4<z<0.75 galaxies observed 
p — , with FLAMES/GIRAFFE at the VLT. Their kinematics had been classified as rotating disks (11 galaxies), perturbed rotators (7 galaxies), or 
f"*) , complex kinematics (14 galaxies). 

Methods. We have computed the specific angular momentum of disks (Jatsk) and the dynamical support of rotating disks through the V/cr ratio. 
To study how angular momentum can be acquired dynamically, we have compared the properties of distant and local galaxies, as a function of 
their kinematical class. 

O ' Results. We find that distant rotating disks have essentially the same properties (jdtsk and Rd) as local disks, while distant galaxies with more 
complex kinematics have a significantly higher scatter in the jdiskr^max. and Rd-Vmax planes. On average, distant galaxies show lower values of 
CZ) ■ V/cr than local galaxies, the lowest V/cr values being reached by distant galaxies showing perturbed rotation. This can probably be attributed to 

' heating mechanisms at work in distant disks. 
J> , Conclusions. We found observational evidence for a non-linear random walk evolution of the angular momentum in galaxies during the last 8 
' ^ n ^ ne ev °l u ti° n related to galaxies with complex kinematics can be attributed to mergers, but not to (smooth) gas accretion alone. If galaxies 
K> ' observed at intermediate redshift are related to present-day spirals, then our results fit quite well with the "spiral rebuilding" scenario proposed 
^ . by Hammer et al. (2005). 

Key words. Galaxies: evolution; Galaxies: formation; Galaxies: kinematics and dynamics; Angular momentum; 3D spectroscopy. 
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1. Introduction 

The evolution and the origin of the Tully-Fisher (T-F) rela- 
tion is still a matter of intense debate. A strong evolution 
both in slope and dispersion has bee n found in B band (e.g., 
Ziegler et al]l2002r. iBohm et al| 2004. from z~l to z=0. More 



recently, IConselice et alJ (120051) derived the T-F relation in K 
band, K-band absolute magnitude being a better tracer of stellar 
mass, and found such a large dispersion that one could wonder 



Flores et al. 



whether or not the T-F relation even exists at z~l 
(2006, tereafter Paper I) derived the first T-F relation for dis 
tant galaxies using integral field spectroscopy. The two dimen- 
sional spatial coverage allowed them to properly identify the 
dynamical nature of distant galaxies, and establish that as much 
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as ~40% of field galaxies are not in equilibrium, i.e., galaxies 
that are not suitable to establish a proper T-F relation. They 
then derived a T-F relation that does not appear to have evolved 
since z~0.6, in the M% (or stellar mass) versus V max plane. On 
the other hand, the stellar m ass density increases from z=l to 
z=0 (e.g. jDrorv et al. I l2005l) : assuming a rough 30% increase 
in stellar mass from z=0.6 to z=0, we should then see a ~0.1 
dex shift along the M» axis between the local and distant T-F 
relations. However, one should keep in mind that typical un- 
certainties usually associated with M* are of the same order, 
i.e., 0.1-0.2 dex (see Paper I). Another possibility is to assume 
some gas accretion from the intergalactic medium which would 
be directly converted into stellar mass: while M* increases, 
Vmax could then increase in the same time, because the total 
mass of the system mainly depends on the rotational veloc- 
ity. Then, both quantities could evolve such as the resulting 
evolution in the T-F plane operat es along the relation (e.g., 



Portinari & Som mer-Larsen 2006). More statistics and/or fu- 
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ture studies of the T-F relation at higher redshifts (where the 
stellar mass density was much smaller than at z=0.6) should 
provide a decisive test of this point. Whether the T-F relation 
evolves with time will no doubt provide important clues about 
relationship between the growth of mass and the characteristics 
of the stellar populations, and, ultimately, about galaxy forma- 
tion and evolution in general. 

The ideas concerning the origin of the T-F relation can 
be divided into two broad categories. In the first one, the T- 
F relation originates from the cosmological equiva lence be- 
tween the halo mass and the circular velocity (e.g., IMo et al.l 
19981) . The relation then comes from the fact that the finite age 
of the Universe imposes a maximal radius from where mat- 
ter can be accreted to form a galaxy. The second broad cat- 
egories of models invoked to explain the T-F relation is self- 
regulated star formation in disks of different mass (e.g., Silk 
1997b . However, numerical simulations taking into account 
both ingredients of gas accretion and self-regulation have not 
been able to reproduce all aspects of the T-F rel ation, such 
as the zero point (e.g.. lSteinmetz & Navarro I ll999h . Many au- 
thors suggested that feedback from massive star formation or 
active gal actic nuclei could help to solve these discrepancies 
(see, e.g., Eisenstein & Loeb 1996; Heavens & Jimenez 1999; 
Steinmetz & Navarro II 19991) . 

Feedback has also been suggested to solve the so-called 
"angular momentum catastrophe" of the A-CDM model, i.e., 
the inability of simulations to reproduce disk galaxies with 
sufficient angul ar momentum in compariso n with what is 
observed (e.g., ISteinmetz & Navarro I Il999h . Feedback has 
been proposed as a potential solution to this problem, which 
has been investigated in detail (e.g., Mailer & Dekel 2002; 

~ 20061) . State-of-the-art 



D'Onghia et al. 1120061 iGovernato et al 



numerical simulations, including the effects of AGN feedback 
driven through accretion of gas onto a super-massive central 
black hole, show how a disk can re-form after the merging 
of two rotating gas-rich disks, with a sufficient and c onsistent 
amount of angular momentum ( Ro bertson et al. 1 l2005l) . 

The mechanisms through which galaxies may have ac- 
quired their a ngular momentum has been discussed for man y 
decades (e.g., IStromberglll934t iHoyle 1 1 195 ft iMestell 1 1963h . 
The so-called "gravitational instability paradigm", indepen- 
dent of the details of the cosmology, predicts that the angu- 
lar momentum of a protogalaxy should grow linearly with 
time due to tidal torques from interactions with neighb oring 



struct u res, until it de couples from the Hubble flow (IPeebles 
White] 119841) . In a more modern picture, galaxies 



form from infalling baryonic gas embedded into dark mat- 
ter haloes, and their angular momentum is then inherited 
from the halo JWhite & R eesl 119781 IF all & Efstathioul ll980t 
Barnes & Efstafhiouf 1987 ). After the end of the epoch dur- 
ing which tidal torquing is effective, subsequent evolution of 
the angular momentum takes place non-linearly, through a 
random walk pr ocess associated with mergers events and/or 
mass accretion JVitvitska et alJ 120021: iMaller & Dekell [2002; 



Peirani et alJl2004l) . This random walk leads to a change in 
the angular momentum of the haloes, with a more signifi- 
cant change (increase or decrease, depe nding on the geome- 
try of the merger) during major mergers dVitvitska et al.l l2002; 



Peirani et al.l 12004). For a given galaxy, it is in this way that 
major mergers are the main source of either positive or nega- 
tive change in its angular momentum. 

Relating the angular momenta of the halo and the disk is 
not at all straightforward. A reasonable assumption often made 
is that the specific angular momentum j (i.e., the angular mo- 
mentum per unit mass), is conserved during the collapse of 
the gas dMestellll963l IFall & Efstathioulll980l) . If we also as- 
sume that bot h gas and dark matte r are well mixed in the 
proto-galaxy jFall & Efstafhi ou 1980), then this leads to jdisk ~ 
jhaio- Making this assumption allows models and simulations 
to reproduce several properties of local disk galaxies (e.. 



Mo et all 1 19981: ISver et all 1 19991: Ivan den Bosch] l20oil 12002 



Toni ni et a 1. 2006). However, some problems remain, such as 



the "mismatch of angula r momentum profiles" between the 



dark matter and the disk (Bullock et al. 2001; van den Bosch 



2002^lMaller & Dekell2002l) . and the inability of models to si 



multaneously match characteristics such as the slope and zero- 
point of the T-F relation, slope and zero-point of the radius- 
luminosity relation, the luminosity function of spirals with rea- 
sonable values for the disk masses, halo structural parame- 
ters, and circular v elocity relative to virial velocity (see, e.g., 
Dutton et alJl2006l) . 



The build-up of angular momentum in rotating disks could 
be better understood through the comparison of theoretical 
models and simulations with observations of distant galax- 
ies. In a first at tempt at estimating the angular momenta of 
distant galaxies. iForster Schreiber et al. (2006) studied several 
z~2 galaxies, and found that, ~10 Gyr ago, galaxies appear to 
have approximately the same specific angular momentum as 
today's spirals, with values roughly similar to that expected for 
their haloes (see also lNesvadba et al. 2006). They argued that 
this confirms the hypothesis that baryons likely acquired their 
angular momentum during the collapse of their parent dark 
matter halos. However, as they pointed out, the spatial reso- 
lution of their observations does not allow to uniquely distin- 
guish between rotating disks and merger-induced kinematics, 
thus the origin of high angular momentum in z~2 disk galaxies 
still remains unclear. 

At lower redshift, a sample of z~0.6 galaxies has 
been observed using the multi-integra l field spect r ograph 



FLAMES/G I RAFF E at VLT (see iFlores et al. I £006; 



Puech et al. 120061 hereafter Paper I and Paper II, respec- 
tively). GIRAFFE observations are confronted with similar 
difficulties in term of spatial resolution as any other study of 
high redshift galaxy dynamics. However, to mitigate against 
these effects, we developed a classification scheme based on 
the kinematics and morphologies of the galaxies, separating 
them into rotating disks and galaxies with complex or dis- 
turbed kinematics and morphologies (see Paper I for details). 
Given all of the oth er comparisons made in subsequent papers, 
e.g., Paper II and iPuech et al. (2006b, hereafter Paper III), 
this classification method appears very robust (see detailed 
discussions in both Paper I and Paper II). The goal of this paper 
is to derive the specific angular momentum in these z~0.6 
galaxies. The kinematics and emission line properties of this 
sample has been studied in detail in the three previous papers 
of this series (dedicated to the GIRAFFE Guaranteed Time 
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Observation sample; see Paper I, Paper II, and Paper III). This 
paper is the forth of this series and is organized as follows: §2 
summarizes the observations and introduce the methodology. 
§3 presents the specific angular momentum of the GIRAFFE 
sample. §4 discusses the dynamical support of distant rotating 
disks, while §5 and 6 discuss the implications and summarize 
our results. 

2. Observations & Methodology 

As part of the Guaranteed Time Observations (GTO) of the 
Paris Observatory, we obtained observations with the multi- 
integral field spectrograph FLAMES/GIRAFFE, of a complete 
sample of 32 galaxies, with redshifts ranging from 0.4 to 0.7, 
EW (OII)> 15A and 1^ <22.5. Briefly, we used the LR04 
and LR05 grating, changing the setups to specifically target the 
[OII]Xl3726,3729 doublet (with R~ 10000); integration times 
ranged from 8 to 13 hours; the seeing was typically ~ 0.6 - 0.7 
arcsec during all the observations. D ata cubes were red uced 
using the GIRBLDRS vl.12 package felecha et al. Il2000l) . in- 
cluding a flat-fielding. Sky was carefully subtracted with our 
own IDL procedures. We derived for these galaxies both veloc- 
ity fields and velocity dispersion maps in Paper I and Paper II. 
These maps, as well as HST morphology, were used to divide 
the sample into three distinct classes based mainly on their dy- 
namical characteristics: rotating disks, perturbed rotators, and 
galaxies with complex kinematics. A complete description of 
the GTO sample and of the methods and analysis used to clas- 
sify the galaxies from our sample are given in Paper I. 

The specific angular momentum j of a rotating system can 
be estimated as (see, e.g. jMo et alll998l) : 



on Rd were directly taken as 1-sigma uncertainties returned by 
GIM2D, with a median value of 0.12 kpc. 



j=J3RV max , 

where B is a dimensionless parameter that depends on the ge- 
ometry and spatial distribution of mass, R a characteristic ra- 
dius of the mass distribution, and W max is the maximal rotation 
velocity. For a thin exponential disk of scale length Rd, this re- 
lation becomes: 

jdisk — ^RdV max- 

To estimate a disk scale length, Rd, it is necessary to de- 
convolve the disk from the bulge componen t. We used GIM2D 
dSimard & Pritchetlll998l:ISimard et alj|2002h to measure R d in 
the galaxies observed with GIRAFFE using HST/WFPC2 im- 
ages (0.1 arcsec pixel -1 ) or with preference, ACS images (0.05 
arcsec pixeF 1 ) when they were available. In Paper I and Paper 
II, we derived the ha lf light radius from th e modeling of isopho- 
tal ellipses (see also lHammer et al.ll2.00lh . Figure Q] shows that 
both half light radii derived using this method and GIM2D 
agree very well. Unfortunately, GIM2D does not estimate any 
uncertainty for the half light radius, but it is noteworthy that 
a linear fit, between the half light radii derived using GIM2D 
and the modelling of isophots, returns a median standard de- 
viation of ~0.29 kpc, in close agreement with the typical un- 
certainty of ~0.34 kpc on the half light radius derived from the 
isophotal ellipses modelling, as claimed in Paper I. Given the 
fitting methods are independent, this comparison suggests that 
our Rrf measurements using GIM2D are robust. Uncertainties 



■Y 




ellipses 



Fig. 1. Comparison of half-light radii (in kpc) derived using 
GIM2D and using the modeled results from fitting isophotal 
ellipses to the galaxy light profile. The resulting half light radii 
derived from both methods agree very well, with the galax- 
ies with disturbed or complex kinematics showing additional 
scatter. Blue dots represent galaxies classified as rotating disks, 
green squares represent perturbed rotations, and red triangles 
represent galaxies with a complex kinematics. 

Unfortunately, for five galaxies, we only had ground-based 
images available w hich were taken at the CFHT (see, e.g., 
Schad e et al. 1996, ; 0.207 arcsec pixel '). Because these im- 
ages are significantly affected by the relatively large seeing disk 
(large compared to the Rd of the disks), we did not attempt to 
use GIM2D to fit their disk light profiles since certainly the 
solution would be highly degenerate, and thus unlikely to be 
reliable. Instead, we used the fact that for an exponential disk, 
Rd can theoretically be obtained from the half light radius rhaif 
using /?d=rhaif/l-68. In Figure|2] we compared Rd as measured 
by GIM2D with niatf/1.68 as deduced from Paper I and Paper 
II. This figure shows that this is a reasonable assumption (at 
least for rotating quiescent disk galaxies). We also used this 
method for two galaxies for which GIM2D failed to provide a 
statistically robust fit to the light profiles. 

Both Figures[T]and[2]show that for galaxies exhibiting more 
complex kinematics than a simple rotating disk (i.e., those clas- 
sified as perturbed rotation or complex), the estimates of Rd 
and/or fhaif nave a larger relative dispersion. In general, how- 
ever, even the galaxies with complex kinematics have Rd and/or 
Hiaif consistent with the average results of the rotating disks. 
The relatively high dispersion can be easily understood as a re- 
sult of their unrelaxed dynamical state: it is likely that these 
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Rhalf ... /1.68 
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Fig. 2. Comparison of disk scale lengths (in kpc) derived using 
GIM2D and using the modeled results from fitting isophotal el- 
lipses to the galaxy light profile. The resulting disk scale length 
from both methods agree very well, with the galaxies with 
disturbed or complex kinematics showing additional scatter. 
Blue dots represent galaxies classified as rotating disks, green 
squares represent perturbed rotations, and red triangles repre- 
sent galaxies with a complex kinematics. The median uncer- 
tainties, i.e., 0.12 kpc for GIM2D, and 0.34/1.68 for the mod- 
elling of the isophots respectively, are indicated in the upper- 
left corner. 



galaxies are not in rotational equilibrium, and an exponential 
disk model is then probably an inadequate representation of 
their true light profile. However, to get homogeneous and con- 
sistent estimates, we choose to treat the whole sample, what- 
ever the kinematical class is, as if all galaxies were rotating 
disks. We also note that despite the increased scatter between 
the various estimates of the scale length and half-light radius, 
each kinematic subsample generally falls along the one-to-one 
line. It is simply that the individual estimates are less reliable 
for the galaxies with complex kinematics but not for the en- 
semble of each class. 

Vmax was estimated by fitting a double gaussian to the [Oil] 
doublet (see Paper I). This introduces a relative random uncer- 
tainty on the measurement of V max that is estimated to be ~ 
10 km/s, from the comparison of independent fits to the data. 
This uncertainty is relatively low because of the high spectral 
resolution of GIRAFFE (R~ 10000) which allows a good ve- 
locity measurement. Moreover, the limited size of the IFU can 
introduce an observational source of uncertainty, if this size is 
too small to reach the flat part of rotations curves. However we 
showed in Paper I that the size of the GIRAFFE IFU is well- 
suited to measure V max for all the galaxies of the sample. 



Because of the coarse sampling of the GIRAFFE IFU data 
(pixel size of 0.52 arcsec, ~ 3.5 kpc at z~0.6), a mean cor- 
rection of 20% has to be applied on V ma x . This correction fac- 
tor is appropriate for z~0.6 rotating disks (see discussion of 
this in Paper I). This correction factor was determined by sim- 
ulating GIRAFFE observations from hydrodynamical simula- 
tions of a Sbc-like rotating disk. As explained in Paper I and 
Paper II, we are probably severely underestimating this correc- 
tion factor for objects showing more complex kinematics (see 
Paper II), and this could introduce a relative offset between ro- 
tating disks and objects with more complex kinematics. The 
absolute uncertainty of the correction factor for rotating disks 
was determined to be +4% because of possible variations in 
size and inclination (see Paper I). An additional source of un- 
certainty on this correction factor could be due to variations 
of the dynamical properties along, e.g., the Hubble sequence. 
We checked that, using another hydrodynamical simulation of 
a "mean" SDSS galaxy, a mean 20% correction factor is still 
appropriate, within + ~2% relatively to 20% correction factor. 
We emphasize that this uncertainty is very difficult to estimate, 
given the high number of parameters which have to be taken 
into account. 

Finally, uncertainties on jdisk were estimated using usual 
methods, i.e. propagating the individual uncertainties on V max , 
Rd, and on the correction factor on V max , as detailed above. 
This gives a median 1-sigma uncertainty of 0.09 dex and 0.10 
dex in V max and jdisk, respectively. This is to be compared with 
typical uncertainties for local galaxies of 0.03 and 0.06 dex re- 



spectively, as estimated from lCourteaul (119971) . We emphasize 
that such an uncertainty is meaningful only for rotating disks, 
since the accuracy and appropriateness of the model used to 
reproduce the data is not taken into account. This uncertainty 
should then be viewed as a random uncertainty associated with 
jdisk, if the rotating disk model is correct, but does certainly not 
include the systematic uncertainty when such a model is not 
appropriate. Indeed, in Papers I and II, we emphasized the fact 
that objects showing complex kinematics are probably merg- 
ers or merger remnants. For these objects, the jdisk derived as 
above should then be viewed as the orbital angular momentum 
due to the relative motion of the two progenitors, rather than 
the intrin sic spin angular momentum of a single rotating disk 
(see also lForster Schreiber et alJl2006l) . We will discuss this in 
more detail subsequently. 



3. Angular Momentum of Distant Galaxies 

To compare distant wit h loca l gal axies, we used two loca l 

d!992h . 



samples from Courteaul ( 1 19971) and iMathewson et al 



For the first sample , we u sed the compilation made by 



Steinmetz & Na varro dl999h . and kindly provided by M. 



Steinmetz. For the second sample, we retrieved the electron- 
ically available data from the CDs3- For this sample, no Rj 
measurement were directly available: we derived it from their 
(I band) 23.5 magnitud e isophotal radius dividi ng by a mean 
correction factor of 3.5 ( Palu nas & 
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Figure [3] shows the disk scale length versus the maximal 
rotation velocity, both for local and distant samples. Both local 
samples agree quite well, in spite of the different proxies used 
for the maximal rotational velocity and the disk scale length. 
Distant rotating disks fall close to local rotating disks, although 
some of them appear to have a relatively lower disk scale 
length. Those are associated to galaxies with half light radius 
slightly lower than that of compact galaxies (i.e., Rhaif < 4.7 
kpc, see Paper II). On the other hand, more kinematically dis- 
turbed distant galaxies show a very high dispersion around the 
local Rd — V max relation. It can be explained in two ways. First, 
more kinematically complex distant galaxies have a larger un- 
certainty in the determination of their radius. This uncertainty 
can be estimated from the scatter of kinematically complex 
galaxies in Figure |2l and is <x ~1.45 kpc. Relatively to the me- 
dian Rd of kinematically complex galaxies, this translates into 
a 0.2 dex scatter in Figure [3] It is thus clear that this is not 
sufficient to explain the extremely large scatter of the kinemat- 
ically complex galaxies in the Rd - V max plane (i.e., cr ~0A5 
dex relatively to the fit of the local galaxies Rd - V max rela- 
tion). The second (and most important) source of scatter is that 
rotational velocities undergo abrupt spatial and temporal vari- 
ations during dynamical processes such as minor and/or ma- 
jor mergers. This effect is illustrated by the black pentagons 
track: they represent a hy drodynami c al sim ulation of a merger 
of two Sbc galaxies from lCox etall d2006l) . GIRAFFE obser- 
vations were simulated along the merging sequence, as well as 
HST images (Puech et al., in preparation). From these simu- 
lated data, both rotation velocities and disk scale lengths were 
extracted and corrected following similar methods used for real 
distant galaxies (see previous section). During this simulation, 
the merger occurs at ~ 1.8 Gyr after the beginning of the simu- 
lation, and is simulated for up to 0.5 Gyrs after the merger. At 
the end of the simulation, the remnant looks relatively similar 
to an elliptical, with a small inner disk at the center. It is be- 
yond the scope of this paper to explore for every possible track 
in Figure[3] varying, e.g., the gas fraction, orbital parameters, 
or the amount of feedback. What is important to note here is 
that this simulation, which includes possible observational or 
instrumental biases, illustrates how the scatter of the most kine- 
matically disturbed distant galaxies could be easily reproduced 
by such events. This strengthens the idea that these objects are 
probably ongoing merger and/or merger remnants. 

Figure H shows the specific angular momentum of the ro- 
tating disks in the sample of galaxies observed with GIRAFFE, 
jdisk, versus their maximum rotation velocity, W max . In this fig- 
ure, z~0.6 rotating disks show a specific angular momentum 
comparable to that of local galaxies with same V max (i.e., ~ 
same total mass). Objects with a complex kinematics have a 
lower median specific angular momentum, with log jdisk ~ 2.7, 
than in rotating disks, which have a median of log jdisk ~ 
3.0. This is probably again due to the fact that V max can be 
severely underestimated for galaxies with complex kinematics 
(see Paper I and Paper II, and above): their derived specific an- 
gular momentum should then be considered as lower bounds. 
As in Figure [3] the dispersion of distant galaxies with com- 
plex kinematics is significantly higher. As these objects are hy- 
pothesized to be undergoing a merger or are merger remnants, 
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Fig. 3. Disk scale length versus maximal rotational velocity, 
for the sample of 32 distant galaxies presented in Paper I 
and Paper II. Blue dots represent galaxies classified as rotat- 
ing disks, green squares those classified as perturbed rotators, 
and red triangles galaxies with complex kinematics . The s mall 
black dots represent the local sample of ICourteaul ( 1 1997b and 
Mathewson et alJ (11992). The median 1-sigma uncertainties are 
indicated in the upper-left corner, for both distant rotating disks 
(blue dots), and local disks (black dots). For local galaxies, Rd 
have been determined using ID fits whereas 2D fits were used 
for distant galaxies: this explains why the uncertainty along the 
ordinate is larger for local than for distant galaxies. Black pen- 
tagons represent simulated GIRAFFE observations using a hy- 
drodynamical simulation of a major merger of two Sbc galax- 
ies. The beginning and the end of the sequence, as well as the 
merger itself, are indicated by solid symbols. 



we are likely measuring an orbital angular momentum rather 
than an intrinsic spin angular momentum of a single disk. In 
such a case, the shape factor B-2 is probably not appropriate, 
and could be an source of additional dispersion. To investigate 
this effect, we also over-plotted in this Figure the same simu- 
lation of GIRAFFE observations of a merger of Sbc galaxies, 
as described above. Note that the beginning of the sequence 
corresponds to only one of the two progenitors, which shows a 
deficit in specific angular momentum compared to local galax- 
ies. This could be related to the difficulty for hydrodynami- 
cal simulations to reproduce the higher angular momentum ob- 
served in local disks. The position of most of galaxies with 
complex kinematics in the jdisk - V max plane can here again be 
easily explained as a result of mergers (see black pentagons). 
As these simulated data include many possible sources of un- 
certainty, it clearly suggests that the larger scatter of galaxies 
showing complex kinematics can be associated with mergers 
or merger remnants. 
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Fig. 4. Specific angular momentum of the disk, jdisk> versus 
their maximum rotation velocity, V max , for the sample of 32 
distant galaxies presented in Paper I and Paper II. Blue dots 
represent galaxies classified as rotating disks, green squares 
those classified as perturbed rotators, and red triangles galaxies 
with complex kinematics . The smal l black dots repre s ent th e 
local sample of lCourteaul (1 1997b and lMathewson et alJ (119921) . 
The median 1-sigma uncertainties are indicated in the upper- 
left corner, for both distant rotating disks (blue dots), and lo- 
cal disks (black dots). Black pentagons represent simulated 
GIRAFFE observations using a hydrodynamical simulation of 
a major merger of two Sbc galaxies. The beginning and the end 
of the sequence, as well as the merger itself, are indicated by 
solid symbols. 

4. Dynamical Support of Distant Disk Galaxies 

In this section, we compare the dynamical support (rotation vs. 
velocity dispersion) of distant and local rotating disks. Such an 
analysis could help us to understand how rotating disks acquire 
and lose their angular momentum. 

The most revealing quantity to estimate the dynamical sup- 
port of galaxies is the ratio of rotation velocity to velocity dis- 
persion of the disk, V/cr. V is a circular velocity (quantifying 
the amount of rotation), and cr is an estimate of the intrinsic 
velocity dispersion in the disk (i.e., turbulent or peculiar mo- 
tions). In elliptical galaxies, V/cr is usually estimated via the 
ratio of the maximal rotational velocity to the mean veloc- 
ity dispersion with in 0.5^, where r f is the effective radi us of 



the system (e.g.. iDavies et al.l 1 19 83: Bend er et al.l 1 1994). For 



spiral galaxies, there is no general consensus: different spa 
tial components of cr estimated fol l owing different methods 
have been used (e. g., Bottemal 1993 ; Vega Beltra n et al. 2001 



sic V/cr is more robustly estimated using the ratio between the 
mean squared rotational velocity and the mean squared velocity 
dispersion, both measured directly along the line-of-sight. 

We choose to estimate V using the maximal rotational ve- 
locity (corrected from inclination and spatial resolution effect, 
see § 2), as it is probably the most accurate quantity derivable 
from GIRAFFE velocity fields. Because the velocity gradient 
of the rotation curve of z ~0.6 rotating disk galaxies falls ap- 
proximatively in only one GIRAFFE IFU pixel (or about one 
resolution element given the seeing, see §2), the center of the 
GIRAFFE velocity dispersion maps show a peak that is due 
to shear and/or large-scale motions in velocity, and cannot be 
used to estimate the intrinsic velocity dispersion of the disk (see 
Paper I and Paper II). On the contrary, the regions surrounding 
the peak of the velocity dispersion are much less affected by 
the shear and/or large-scale motions in velocity, since the ro- 
tation curve is approximatively flat in these regions (i.e., has a 
constant velocity with radius). These outer regions of the ve- 
locity dispersion maps can then be used to construct a reliable 
estimate of the intrinsic velocity dispersion of the disk. To con- 
struct such reliable estimates of the intrinsic velocity dispersion 
of the disks, we first removed the cr peak due to rotation in the 
velocity dispersion map, guided by the modeling of the disper- 
sion we made in Paper I. We then estimated cr by deriving the 
signal-to-noise weighted mean of the remaining pixels. 

Unfortunately, there are very few published velocity dis- 
persion maps of local galaxies that could be used for direct 
comparison. However, several dozen velocity dispersion pro- 
files of the gas in local disks that were obtained with long- 
slit spect roscopy have been publis h ed. We combi n ed the sam- 
ples of IVega Beltran etall (boOll) , ICorsini et ail d2003l) . and 
Pizzel la et al.l d2004 . which are composed of spiral galaxies 
with morphological type earlier than Sc. Note that V/cr ratio 
does not seem to depend on the morphological type, since a 
similar range of values are foun d for both spiral and irregular 
galaxies (see iHunter et alj|2005l) . We checked that both local 
and distant samples have similar distribution in the absolute 
B magnitude, Mb, and we kept only galaxies with absolute B 
magnitude brighter than the lower value found in the GIRAFFE 
sample (i.e., Mb(AB) < -19.26), to ensure that both samples 
are probing galaxies with comparable stellar masses and star- 
formation rates. To exclude the central dynamically hot region, 
we performed a sigma clipping on the velocity dispersion pro- 
file of local galaxies, keeping only the points along the curve 
that were below 2cr around the median. We then took the final 
median of the remaining points as a measure of cr in the disk. 

When measured projected onto inclined disks, the observed 
velocity dispersion is a combination of the three spatial compo- 
nents of the velocity dispersion - the radial component cr,., the 
azimuthal component cr^, a nd the vertical component cr. (e.g., 
Binney & Merrifield Ill989l) : 



Hunt er et al.ll2005l) . iBinneyl (120051) demonstrated that with in- 
tegral field spectroscopy (i.e., a 2-dimensional spatial cover- 
age combined with simultaneous spectral coverage), the intrin- 



cr 2 ohs - (crj sin 2 r] + cr 2 cos 2 77) sin 2 i + cr\ cos 2 i 

where 77 is the angle between the observed PA and the ma- 
jor axis of the galaxy, and i is the inclination angle of the 
disk. In local spira l galaxies, cr z ~ cr^ ~ 0.7cr r for stars (e.g., 
Hunter et al] |2005). In the case of local galaxies, if we assume 



both a well aligned slit and that gas and stars dynamics are well 
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coupled (see lVega Beltran et al. 2001; Pizze lla et al .12004 . one 
gets cr ohs ~ cr,. In the case of distant galaxies, observed us- 
ing integral field spectroscopy, it is necessary to correct for 
the different 77 corresponding to the different IFU pixels. We 
can directly correct for this effect on the mean sigma by av- 
eraging the cos 2 1] and sin 2 77 terms, and then multiplying by a 
(1 + sin 2 2'/2) _1/ ' 2 correction factor. Finally, we obtain an es- 
timate of the spatially averaged cr z , that can be used in the 
V/cr ratio. To check that the final estimate of V/cr does not de- 
pend significantly on the method used to estimate cr, we com- 
pare d with the velocit y dispersion measured at Hiaif/4 as given 
by IVega Beltran et al! . Using this alternative method does not 
significantly affect the general trend of V/cr estimated in local 
disks. 

In the sample of local galaxies, errors have been estimated 
as follows. The uncertainty on V was directly taken as the mean 
of the ve locity measuremen t uncer t ainty on V max and V m j n as 



30 



given by [V ega Bel tran et all d200lb , ICorsini et all d2003l) . and 
Pizzella et al.1 d2004l) . The uncertainty on cr was estimated as 
the median of measurement uncertainties of the points of the 
velocity dispersion curve kept during the sigma clipping (see 
above). A conservative uncertainty of three degrees was as- 
sumed for the inclination. We found a median uncertainty of 
1.6. Figure lATI shows how this uncertainty evolves with V/cr. 
This Figure shows that high values of V/cr have large associ- 
ated uncertainties; including only sources with V/cr estimates 
greater than three times their uncertainty, left us with a range 
of local V/cr values ranging between 2 and 9, with a median 
uncertainty in V/cr of 1.1. In the following, we will consider 
these points only, as others galaxies have V/cr ratios that are 
too uncertain. For distant rotating disks, we estimated the un- 
certainty of cr to be lower than 15 km/s (see Paper I). We added 
an additional uncertainty term due to the correction factor used 
to deconvolve V max (see § 2). 

Figure [6] shows the comparison of V/cr as a function of the 
ellipticity in both distant and local rotating disk galaxies. The 
median V/cr is 3.8+2 for distant rotating disks, and 6.1±1.1 
for local disks. Recall that highest local V/cr have been dis- 
carded in this comparison because of a too high uncertainty: 
this could result in an underestimation of the local median 
V/cr. We also plotted the V/cr ratio for the distant perturbed 
rotating disks, with a median of 2.4+2.5. Figure [6] shows that 
distant disks have lower V/cr ratios than local ones. This is 
consistent with the fact these galaxies are likely undergoing 
a minor merger and/or a gas accretion event which is heat- 
ing their disks (W alker et al. 1996; IVelazquez & Whitd 1999; 
Abadietal 2003), as was suggested in Paper II. More data will 
be needed to confirm this trend, however. 

5. Discussion 

5.1. The Growth of Angular Momentum in Galaxies 
with Time 

We showed in § 3 that z~0.6 and local disk of same total mass 
show comparable specific angular momentum. More dynami- 
cally complex objects fall around the same jssk - V max relation 
but with much larger dispersion than the local relation. How 
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Fig. 5. V/cr vs. uncertainty in V/ cr for the sample of local galax- 
ies used for comparison. The galaxies with significant V/cr es- 
timates (V/cr > 3 A (V/cr)) are shown in red. 



can we interpret this in terms of angular momentum growth 
with time in galaxies ? 

One of the first theories proposed for the build up of angu 



lar momentum in g alaxies is the tidal torque theory dPeebles 



19691 IWhite II 19841) . The basic idea in the tidal torque theory 
is that most of the angular momentum is being gradually ac- 
quired in the linear regime of growth of the density fluctua- 
tions, due to tidal torques from neighboring fluctuations. This 
process continues relatively efficiently until the halo reaches 
its turn around time, i.e., roughly when the protogalaxy de- 
couples from the Hubble flow. Moreover, this theory suggests 
that angular momentum gain is minimal when haloes are grow- 
ing non-linearly after they have decoupled from the expand- 
ing background and have formed virialized systems. This the- 
ory then predicts rapid growth of angular momentum in the 
early evolution of m assive halos, and subsequently little growth 
( Porciani et alJl2002|). At first sig h t, our result s, combined with 
those of iForster Schreiber et al. (2006) and iNesvadba et al. 



( 2006) would mean that the specific angular momentum of ro- 
tating disks does not evolve with redshift from z~2 to present. 
This would indicate that this scenario for the growth of an- 
gular momentum is plausible. However, we have to caution 
that current samples are relatively small and incomplete in 
terms of galaxy mass. Moreover, there are some issues with 
the physical resolution of the data, especially at z>2 (but see 
Nesvadba et al. 1 120060 . Even though these results are all for- 
mally consistent with no overall evolution in the specific angu- 
lar momentum, given the limitations of all studies, it is difficult 
to conclude anything about evolutionary trends in the specific 
angular momentum in rotating disk galaxies. 
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Fig. 6. Comparison of V/cr and ellipticity, e, for galaxies at 
z~0.6 compared to a matched sample of local galaxies. The 
V max /o~ ratio for galaxies at z~ 0.6 that have been determined 
to be rotating disks (blue dots) and to have rotation curves that 
are perturbed (green squares) are compared with local galaxies 
whose morphology and velocity curves are consistent with ro- 
tating disks (black crosses). The medians are indicated by the 
position of the median uncertainties on the right side of the di- 
agram (each being represented by matching colors) 



Beyond this, the results presented here can constrain 
how angular momentum and indeed how galaxies grow af- 
ter the epoch when massive halos (like the ones studies 
here) acquired their angular momenta through tidal torquing. 
The basic ideas that have emerged from theoretical research 
are that the specific angular momentum can be increased 
or decreased in a halo through a random walk p rocess of 
interactions and merging with other halos (e. g. 



2001; Vitvitska et al 



Gardner 



2002; IPorciani et all l2002h or thro ugh 



orientation of its angular momentum vector relative to the per- 
haps pre-existing disk, etc. 

Returning to both Figs. [3] and [4] it is clear that perturbed 
rotation and kinematically complex galaxies are dynamically 
different from dynamically relaxed local disks. The dispersions 
in the jdisk-Vmax are indeed similar for the rotating disks, higher 
for the perturbed rotators, and very high for the galaxies with 
complex dynamics. We can also note that some of the dis- 
tant galaxies which are not rotating disks have offsets to both 
lower and higher jdisk relative to local galaxies for the same 
V max . Because we are likely underestimating jdisk in kinemati- 
cally complex galaxies (see § 3.1), we cannot conclude specif- 
ically about any difference on median jdisk between the dif- 
ferent kinematical classes. However, the higher dispersion in 
the V max - jdisk plane of kinematically complex galaxies rela- 
tively to distant rotating disks is significant. This dispersion 
is caused by abrupt variations of the angular momentum in 
galaxies with complex kinematics, and are likely the result of 
a random walk during a non-linear phase of evolution of the 
angular momentum. This is the first observational evidence for 
such a non-linear growth of the angular mom entum in galax 



ies, as expected from theoretical m odels (e.g.,|Gardner 
Vitvitska et aIll2002tlPeirani et al.ll2004 . 



2001 



An important question is what is the driving mechanism of 
this random walk. Theoretical models show that major merg- 
ers cause the most abrupt variations of the angular momentum, 
while minor mergers and/ or gas accretion are associated with 



smoo t her variations (e.g.. I Vitvitska et al.1 l2002t iPeirani et aL 
2004); iHetznecker & Burkertll2006l) . If the gas is heated to ap- 



gas accretion onto the halo (e.g., Whit e & Reesl 119781) . In 
the merger process, angular momentum is gained or lost 
depending on the geometry, dynamics, the detailed mass 
distribution of the merging halos, strength of the feed- 
back, the relative masses of the merging halos and galaxies , matics is consistent with the merger scenario. IVitvitska et al 



proximately the virial temperature o f the halo before collapsin g 
to form or onto a pre-existing disk (De kel & Birnboim 2006), 
we would expect the gas to have a specific angular momen- 
tum similar or greater than that of the dark matter halo (e.g., 
Chen et al.ll2003t lOkamoto et al.l l2005). This is why we might 
then expect that the specific angular momentum of the disk to 
remain constant or to mildly increase or decrease at late times 
(Pei rani et al]l2004l) . It is hard to believe that smooth gas ac- 
cretion alone could explain the complex or perturbed kinemat- 
ics and/or discrepant values of the specific angular momentum. 
We instead suggest that mergers (both major and minor merg- 
ers, possibly associated with complex kinematic galaxies, and 
perturbed rotations, respectively, see Paper I and II) play a sig- 
nificant role in changing the angular momentum of galaxies 
with time. While overall consistent with the tidal torque the- 
ory, the increased dispersion in specific angular momentum and 
spin of the perturbed rotators and galaxies with complex kine 



etc. (e.g ., Robertson et al. 2005; Sprin gel & Hernquis t 2005 



D'Onghia et al. 1 12006 ). Alternatively, the angular momentum 



of disk galaxies could grow through the accretion of cooling 
halo gas or smooth accretion of mass into the dark matter halo 
(e.g., IWhite & Reeslll978t IPeirani et al.ll2004 Ivan den Bosch 
2001. |2QQ2al Ivan den Bosch et al.1 l2002bl 120031; IChen et al 



2003t lOkamoto et al.1 120051) . The amount of angular momen- 



tum change then depends on the relative angular momentum of 
the gas and dark matter, how the gas was accreted, how much 
mechanical dissipation the gas undergoes during collapse, the 



( 2002) show that within the context of merging, we would find 
both increases and decreases in the specific angular momen- 
tum, resulting in an increase of the dispersion of their angular 
momenta. This is consistent with what we observe when we 
compare distant rotating disks to galaxies with complex kine- 
matics. This interpretation is strongly supported by a simula- 
tion of GIRAFFE observations of a major merger. This sim- 
ulation illustrates how the higher dispersion of kinematically 
complex galaxies can arise from mergers, taking into account 
all possible observational uncertainties . 
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5.2. The Building of Local Disk Galaxies 

In the following discussion, we will assume that most of the 
observed distant galaxies with emission lines are progenitors 
of local disks. This should be the case even for most galaxies 
showing complex kinematics and perturbed rotations. Indeed, 
those galaxies represent 40% of galaxies at z~0.6, and if they 
were E/SO progenitors, there would be a much higher fraction 
of E/SO that it is observed today (see, e.g., Hammer et al. 2005, 
and also Lotz et al. 2006). 

Because we are comparing samples of galaxies span- 
ning different total mass (V mflA ), linking distant to local 
disks from Figure [4] alone is not straight-forward. Can z~0.6 
gala xies evolve towards local disks through major merg- 
ers/ lHammer et all d2005l) claim that 75% of local spiral in 
the 1Q 105 -10 1l 5 Mq range (the so-ca lled "intermediate-mass" 



range, see also lHammer et al.l 120061) have experienced a ma- 
j or merger since z=l . Assuming an evolution rate as (1+z) 27 
dLe Fevre et al. 20001) . one can derive that ~ 29% of local 



spiral gal axies have e xperie nced a major merger since z=0.6. 
Similarly, Lotz et al. (2006) estimated that between 33% and 
66% of Lg > 0AL* B galaxies had a major merger since z=l.l 
with an evolution rate of (1+z) 112 . Combining the above es- 
timates implies that ~ 15 to 30% of local spirals could have 
experienced a major merger since z=0.6. Then, the majority 
of z~0.6 rotating disks cannot evolve towards local spirals 
through major mergers. Another possibility, is that z~0.6 galax- 
ies could evolved to z=0 spirals through minor mergers and/or 
gas accretion: their rotating disks would then have survived, 
and their specific angular momentum should statistically in- 
crease by only ~ 0.1 dex (iPeirani et al.l 120041) . A third possi- 
bility is that distant disks would evolve in a "closed box". In 
absence of external torques, their specific angular momentum 
would then remain constant. Finally, the last two possible evo- 
lution tracks from z~0.6 to z=0 are equally viable. 

This picture is also supported by the fact that distant 
disks (we consider here both rotating disks and perturbed 
rotating disks together) seem to be heated relatively to lo- 
cal disks (i.e., have systematically higher V/ cr), possibly 



through minor mergers and/or gas accretion (|Walker et al 



1995 IVelazauez& White! Il999t lAbadi et all Eo03h . A diffi 



culty to quantify uncertainty in this result could arise from 
the difference observational strategies used in the local and 
distant samples (slit vs. integral field spectroscopy, respec- 
tively). This, however, is unlikely given the fact that most spi- 
rals today are strongly dom inated by rotational motions (e.g., 



Binnev & Merrifieldl[l989l) . An interesting additional compo- 
nent of the evolution in V max /cr we observed is the possi- 
bility of "angular momentum mi xing" during minor mergers 
and/or gas accretion events (see lOkamoto et al.ll2005l) . During 
such events, the direction of the angular momentum vector can 
change, and newly accreted gas may then settle to a differ- 
ent orientation from the pre-existing disk. During this process, 
the observational manifestations are likely to be relatively high 
velocity dispersion of the gas caused by the gas being shock 
heated by turbulence generated by the overlapping and dynam- 
ically different disks (what we could observe in Figure[6]l. 



Given the above discussion (see § 5.1) about the specific 
angular momentum, it is unlikely that (smooth) gas accretion 
alone plays an important role in the general growth of rotating 
disks. It is also equally unlikely that gas accretion alone could 
explain the disturbed kinematics of perturbed rotation, as this 
process is expected to be dynamically smooth (see Paper I and 
II) and only lead to small changes in the angular momentum. 
Note, however, that, to our knowledge, there is no clear theoret- 
ical prediction about the kinematics of gas accretion by galax- 
ies. Thus, our results show that it is likely that mergers (and, a 
priori, both major and minor mergers) play an important role 
in the general build-up of rotating disks. The link between dy- 
namical processes such as mergers, and the dynamical classi- 
fications adopted in this series of paper will be addressed in a 
forthcoming article. 

The results presented in this paper fit quite nicely within th e 
"spiral rebuilding" scenario proposed by lHammer et al.l (120051) . 
where ~75+25% of local spirals (those of early type) have un- 
dergone a merger since z= 1 and have rebuilt a disk thanks to 
gas accretion. This scenario is composed of 3 major phases: 
a "pre-merger phase" during which two distant spirals merge, 
the "LCG phase" where all material from the progenitors fall 
into the mass barycenter of the system and form a bulge, and 
the "disk growing phase" where subseq uently accreted mate 



rial forms a rotating disk (see Fig. 8 of Ok amoto et al.l 12005 



for how this might look). As the two progenitors are merg- 
ing, their disks, and thus, thei r spin angular m omentum, are 
destroyed during the collision ( Cox et al. 2004). At the same 
time, the encounter has a significant quantity of orbital angu- 
lar momentum which can be progressively converted into spin 
angular momentum as a new disk is re-built around the rem- 
nant. Moreover, we saw that the increase of V/cr in disk galax- 
ies could be due to heating mechanisms such as minor merg- 
ers and gas accretion events, as if some of the disk-rebuilding 
were trig gered by accretion of the gas l eft over from the merg- 
ing event. Yoachi m & Dalcanton ( 20061) suggest that the prop- 
erties of thick and thin disks are consistent with gas rich merg- 
ers playing a significant role in their formation where the stars 
in these mergers formed the thick disk while the settling gas 
formed much of the thin disk. The results presented here am- 
plify these ideas and show their credibility in directly explain- 
ing the dynamics of intermediate redshift disk galaxies. 

6. Conclusions 

We have studied the angular momentum and the dynamical 
support of a sample of z~0.6 galaxies observed with the inte- 
gral field spectrograph FLAMES/GIRAFFE. We found that the 
classification of distant galaxies based on their kinematic prop- 
erties (mainly) and morphologies into three distinct classes, 
i.e., rotating disks, perturbed rotators, and kinematically com- 
plex, is apparently robust. This classification appears to also 
select galaxies with angular momenta consistent with local spi- 
ral galaxies but show varying degrees of dispersion relative to 
the local values, increasing from rotating disks to kinematically 
complex galaxies. 

This can be interpreted as an evidence for a non-linear ran- 
dom walk evolution of the angular momentum during the last 
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8 Gyr. A natural driver for this random walk is provided by 
major mergers, since the dispersion of kinematically complex 
galaxies in the - V max plane, as well as the complexity of 
their velocity fields itself (see Puech et al. 2007, in prepara- 
tion), can both be reproduced by hydrodynamical simulations 
of such events. Major mergers also explains how the angular 
momentum of local rotating disks could be acquired and be as 
high as observed, namely through the conversion of orbital an- 
gular momentum to spin momentum. Moreover, distant disks 
appear to be more turbulent (lower V/cr). This is likely an in- 
dication that local disks could grow through the accretion of 
gas through mergers or in discrete clouds. In a refinement of 
this general picture, the higher velocity dispersions could be 
due to "angular momentum mixing" (see lOkamoto et ai]|2005 ) 
whereby the relative orbit of the infalling gas is skewed com- 
pared to the previous disk. The torque provided by the accreting 
gas causes a change in the angular momentum vector thereby 
increasing the dispersion observed in the gas. 

Finally, these fin dings are consistent w ith the observational 
scenario proposed bv lHammer et al. (2005), where ~15% of lo- 
cal spirals (those of early type) have undergone a major merger 
since z= 1 , and have then rebuilt their disks thanks to gas accre- 
tion (possibly from high angular momentum gas left over from 
the merging event) and/or minor mergers. Unfortunately, draw- 
ing robust conclusions as to the physical processes driving the 
observational manifestations we have discussed in this paper is 
difficult given the small numbers of galaxies we have studied. 
To increase the statistical robustness of these results, we are 
currently analyzing similar data for a much larger sample of 
several hundred disk galaxies as part of the VLT Large Program 
IMAGES (PL: F. Hammer). The galaxies are selected by abso- 
lute J-band magnitude, to have redshifts of or less than 0.9, 
and to have fOUl equivalent width comparable to the galaxies 
study here dRavikumar et al. 112006 ). With the robust selection 
and large numbers, we should be able to make more definitive 
statements about the physical behind the dynamical and mor- 
phological evolution of spiral galaxies over the last 7 Gyrs. 
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Appendix A: Spin Parameter of the Haloes 

We present in this appendix a perhaps more predictive com- 
parison between local and distant galaxies' angular momen- 
tum, based on the use of the spin parameter of the sur- 
rounding haloes, which has the advantag e of being approx- 
imate l y independent of the halo m a ss ([Barnes & Efstathiou 

M 



19871: iLemson & Ka uffmann 1 19991: iMaller & Dekell 1200: 



Peir ani et alJ 12004). This dimension-less parameter measures 



the ratio between the true angular momentum of the system 
and the angular momentum it would have, if all the mass was 
entirely supported by rotation, or, equivalently, the ratio of the 
spin energy to the total binding energy. The spin parameter 
(Peebles 1969) is given by: 



A = 



m 



1/2 



GM 5 ' 2 ' 



where J is the true halo angular momentum, E the total en- 
ergy and M the total mass of the system. It is possible to es- 
timate A using observe d characteristics of the disk. Follow ing 
the method outlined bv lHernandez & Cervantes- Sodi (2006), ,1 
can be estimated using: 

A ~ R d /V 2 d . 

This formulation requires the following assumptions: (1) 
jdisk ~ halo ( see discussion in the introduction); (2) the total 
disk mass must be a constant fraction of the mass of the halo; 
(3) the total energy is dominated by the halo which is virial- 
ized; (4) a baryonic Tully-Fi sher relation e xists which can be 
expressed by M c i ~ Vj (see lMcGaughll2005lH We additionally 
require for our analysis to be valid that these assumptions are 
not dependent on the epoch of observation. In Figure lATl we 
plotted the histogram of A in the local sample used in this pa- 
per. The dash line at A ~ 0.0 3 is the most probable local value 
found by lTonini et al. ( 20061) using an independent method. It 
is very close the peak of our histogram: in spite of its crude- 
ness, this suggests that the method to estimate A is appropriate, 
at least for local disks. 



2 iHernandez & Cervan tes-Sodi ( 2006) used a baryonic Tully-Fisher 
relation slightly different of the one used in this paper, leading them 



toA~R d /V 3 d 5 . 
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Fig. A.l. Histogram of A derived using the method outlined in 
the text for the local sample of gala xies. The most prob able 
value is very close to estimate made by lTonini et al.l (12006) 
using a local sample of galaxies. 



i also 



In Figure DOl we used this relation to investigate the differ- 
ence in th e spin parame t er, A , between local and distant galax- 
ies. From iPeirani et al.l (120041) . the median A between z=0.6 
and z=0 is expected to remain roughly constant with accretion 
only, whereas an increase of ~10% is expected, taking into ac- 
count major merger events. Figure \A2\ seems to favor a small 
increase in A from z~0.6 to z=0. Unfortunately, the uncertainty 
on the median A of the z~0.6 galaxies is very high, see Fig. 
IA.2I it is then not possible to identify the relative contribution 
of these two mechanisms. Both perturbed rotations and objects 
with complex kinematics show higher dispersion in the spin pa- 
rameters than rotating disks (both local and distant). Part of this 
dispersion is probably due to the rotation velocity being under- 
estimated for the galaxies which are not dynamically relaxed 
(see above). It could be tempting to claim for an increasing 
dispersion among the three dynamical classes, from rotating 
disks to objects with complex kinematics, as was seen in the 
specific angular momentum, in agreement with a random walk 
model driven by mergers. However, uncertainty also increases 
from rotating disks to kinematically complex galaxies, and it 
becomes then very difficult to firmly claim for such an evolu- 
tion in the dispersions. Although it seems possible to measure 
A for local galaxies, its derivation seems still very difficult for 
more distant objects. 



Fig.A.2. Spin parameter, A, vs. ellipticity, e, in local and dis- 
tant galaxies samples. Blue dots represent galaxies classified 
as rotating disks, green squares represent those classified as 
perturbed rotators, and red triangles represent galaxies with 
complex kinem atics . Black dots represent th e local sample of 
Courte au ( 1997) and Math ewson et al. (1992). The medians are 



indicated by the position of the median uncertainties for distant 
galaxies are indicated on the right side of the plot (each being 
represented by matching colors). Errors include both standard 
propagation of errors, and a bootstrap estimation of the error 
associated with the determination of the median itself. 



